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Abstract 

We show that the phenomenology of isospin effects on heavy ion reactions at 
intermediate energies (few AGeV range) is extremely rich and can allow a 
"direct" study of the covariant structure of the isovector interaction in a high 
density hadron medium. We work within a relativistic transport frame, beyond 
a cascade picture, consistently derived from effective Lagrangians, where isospin 
effects are accounted for in the mean field and collision terms. We show that 
rather sensitive observables are provided by the pion/kaon production (ir~ /ir + , 
K°/K+ yields). Relevant non-equilibrium effects are stressed. The possibility 
of the transition to a mixed hadron-quark phase, at high baryon and isospin 
density, is finally suggested. Some signatures could come from an expected 
"neutron trapping" effect. 



1 Introduction 



Recently the development of new heavy ion facilities (radioactive beams) has 
driven the interest on the dynamical behaviour of asymmetric matter, H| . Here 
we focus our attention on relativistic heavy ion collisions, that provide a unique 
terrestrial opportunity to probe the in-medium nuclear interaction in high den- 
sity and high momentum regions. An effective Lagrangian approach to the 
hadron interacting system is extended to the isospin degree of freedom: within 
the same frame equilibrium properties (EoS, ^) and transport dynamics can 
be consistently derived. 

Within a covariant picture of the nuclear mean field, for the description 
of the symmetry energy at saturation (04 parameter of the Weizsaccker mass 
formula) (a) only the Lorentz vector p mesonic field, and (b) both, the vector 
p (repulsive) and scalar S (attractive) effective fields ^ E} can b e included. 
In the latter case the competition between scalar and vector fields leads to a 
stiffer symmetry term at high density ^ HI . We present here observable effects 
in the dynamics of heavy ion collisions. We focus our attention on the isospin 
content of meson production. We finally show that in the compression stage of 
isospin asymmetric collisions we can even enter a mixed deconfined phase. 

2 Relativistic Transport 

The starting point is a simple phcnomenological version of the Non-Linear 
(with respect to the iso-scalar, Lorentz scalar a field) effective nucleon-boson 
field theory, the Quantum-Hadro-Dynamics 3 . According to this picture 
the presence of the hadronic medium leads to effective masses and momenta 
M* = M + S s , k*^ = W - W, with S a , E* 4 scalar and vector self-energies. 
For asymmetric matter the self-energies are different for protons and neutrons, 
depending on the isovector meson contributions. We will call the corresponding 
models as NLp and NLp5, respectively, and just NL the case without isovector 
interactions. For the more general NLpd case the self-energies of protons and 
neutrons read: 



E a (p,7l) = -f a v{Ps) ±f5ps3, 



(1) 



(upper signs for neutrons), where p s = p sp + p sn , j a = j% + j%, p s3 = p sp - 
Psn i is = jp — jn are the total and isospin scalar densities and currents and 
fa,u,p,S are the coupling constants of the various mesonic fields. o~(p s ) is the 
solution of the non linear equation for the a field ^ H , 

For the description of heavy ion collisions we solve the covariant trans- 
port equation of the Boltzmann type within the Relativistic Landau Vlasov 
(RLV) method, using phase-space Gaussian test particles and applying 
a Monte-Carlo procedure for the hard hadron collisions. The collision term 
includes elastic and inelastic processes involving the production/absorption of 
the A(1232MeV) and N*(lA40MeV) resonances as well as their decays into 
pion channels, ® . 

3 Isospin effects on pion and kaon production at intermediate en- 
ergies 

Kaon production has been proven to be a reliable observable for the high density 
EoS in the isoscalar sector ^ E| Here we show that the K°' + production (in 
particular the K° /K + yield ratio) can be also used to probe the isovector part 
of the EoS, 

Using our RMF transport approach we analyze pion and kaon production 
in central 197 Au + 197 Au collisions in the 0.8 — 1.8 AGeV beam energy range, 
comparing models giving the same "soft" EoS for symmetric matter and with 
different effective field choices for E sym . We will use three Lagrangians with 
constant nucleon- meson couplings (NL... type, see before) and one with density 
dependent couplings {DDF, see ©), recently suggested for better nucleonic 
properties of neutron stars . 

Fig. [T]reports the temporal evolution of A ±,0 ' ++ resonances, pions i^' ) 
and kaons (K + '°) for central Au+Au collisions at lAGeV. It is clear that, 
while the pion yield freezes out at times of the order of 50/m/c, i.e. at the 
final stage of the reaction (and at low densities), kaon production occur within 
the very early (compression) stage, and the yield saturates at around 20/m/c. 
From Fig. [T]wc sec that the pion results are weakly dependent on the isospin 
part of the nuclear mean field. However, a slight increase (decrease) in the 
7r~ (7r + ) multiplicity is observed when going from the NL (or DDF) to the 
NLp and then to the NLpS model, i.e. increasing the vector contribution f p in 
the isovector channel. This trend is more pronounced for kaons, see the right 
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Figure 1: Time evolution of the A ±,0,++ resonances and pions tt^' (left), 
and kaons (K +, ° (right) for a central (b = fm impact parameter) 
Au+Au collision at 1 AGeV incident energy. Transport calculation using 
the NL, NLp, NLpS and DDF models for the iso-vector part of the 
nuclear EoS are shown. The inset contains the differential K°/K + ratio 
as a function of the kaon emission time. 



panel, due to the high density selection of the source and the proximity to the 
production threshold. Consistently, as shown in the insert, larger effects are 
expected for early emitted kaons, reflecting the early N/Z of the system. 

When isovector fields are included the symmetry potential energy in 
neutron- rich matter is repulsive for neutrons and attractive for protons. In 
a HIC this leads to a fast, pre-equilibrium, emission of neutrons. Such a 
mean field mechanism, often referred to as isospin fractionation is respon- 
sible for a reduction of the neutron to proton ratio during the high density 
phase, with direct consequences on particle production in inelastic NN colli- 
sions. 

Threshold effects represent a more subtle point. The energy conservation 
in a hadron collision in general has to be formulated in terms of the canon- 
ical momenta, i.e. for a reaction 1 + 2 — > 3 + 4 as Si n = (k^ + ht^) 2 = 
(fcg + k±) 2 = s ou t- Since hadrons are propagating with effective (kinetic) 
momenta and masses, an equivalent relation should be formulated starting 
from the effective in-medium quantities k*^ = k^ — S M and m* — m + S s , 
where E s and S M are the scalar and vector self-energies, Eqs.(TT]). The self- 
energy contributions will influence the particle production at the level of thresh- 



olds as well as of the phase space available in the final channel. In fact 
the threshold effect is dominant and consequently the results are nicely sen- 
sitive to the covariant structure of the isovector fields. At each beam energy 
we see an increase of the ir~ /tt + and K°/K + yield ratios with the models 
NL -> DDF -» NLp -> NLp5. The effect is larger for the K°/K+ compared 
to the 7r~ /tt + ratio. This is due to the subthreshold production and to the fact 
that the isospin effect enters twice in the two-step production of kaons, see 
Interestingly the Iso- EoS effect for pions is increasing at lower energies, when 
approaching the production threshold. 

We have to note that in a previous study of kaon production in excited nu- 
clear matter the dependence of the K°/K + yield ratio on the effective isovector 
interaction appears much larger (see Fig. 8 of ref. ©). The point is that in the 
non-equilibrium case of a heavy ion collision the asymmetry of the source where 
kaons are produced is in fact reduced by the n — > p "transformation" , due to 
the favored nn — > pA~ processes. This effect is almost absent at equilibrium 
due to the inverse transitions, see Fig. 3 of ref. ®. Moreover in infinite nuclear 
matter even the fast neutron emission is not present. This result clearly shows 
that chemical equilibrium models can lead to uncorrect results when used for 
transient states of an open system. 

4 Testing Deconfinement at High Isospin Density 

The hadronic matter is expected to undergo a phase transition to a deconfincd 
phase of quarks and gluons at large densities and/or high temperatures. On 
very general grounds, the transition's critical densities are expected to depend 
on the isospin of the system, but no experimental tests of this dependence have 
been performed so far. In order to check the possibility of observing some 
precursor signals of a new physics even in collisions of stable nuclei at inter- 
mediate energies we have performed some event simulations for the collision of 
very heavy, neutron-rich, elements. We have chosen the reaction 238 [/ _|_ 238 JJ 
(average proton fraction Z/A = 0.39) at I AGeV and semicentral impact pa- 
rameter b = 7 fm just to increase the neutron excess in the interacting region. 
In Fig. [2] we report the evolution of momentum distribution and baryon den- 
sity in a space cell located in the cm. of the system. We see that after about 
10 fm/c a local equilibration is achieved. We have a unique Fermi distribution 
and from a simple fit we can evaluate the local temperature (black numbers 



Figure 2: 238 f/ -|_ 238 JJ ^ \ AGeV, semicentral. Correlation between den- 
sity, temperature, momentum thermalization inside a cubic cell, 2.5 fm 
wide, in the center of mass of the system. 

in MeV). We note that a rather exotic nuclear matter is formed in a transient 
time of the order of 10 fm/c, with baryon density around 3 — 4po, temperature 
50 — 60 MeV, energy density 500 MeV fm~ 3 and proton fraction between 0.35 
and 0.40, likely inside the estimated mixed phase region. 

In fact we can study the isospin dependence of the transition densities 
The structure of the mixed phase is obtained by imposing the Gibbs conditions 
for chemical potentials and pressure and by requiring the conservation of 
the total baryon and isospin densities 



where x 1S the fraction of quark matter in the mixed phase. In this way we get 
the binodal surface which gives the phase coexistence region in the (T, pB, pa) 
space. For a fixed value of the conserved charge p$ we will study the boundaries 
of the mixed phase region in the (T,ps) plane. In the hadronic phase the 
charge chemical potential is given by pz = 2E sym (pB) J ^ ■ Thus, we expect 
critical densities rather sensitive to the isovector channel in the hadronic EoS. 




pW(T )A g) = pW)(T, M g), 
Pb = (1 - X)Pb + XPb • 

P3 = (1 - X)P3 + XP3 , 



(2) 
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Figure 3: Variation of the transition density with proton fraction for 
various hadronic EoS parameterizations. Dotted line: GM3 RMF- 
model EJ); dashed line: NLp ; solid line: NLp5 . For the quark 
EoS: MIT bag model with 5 1/,4 =150 MeV. The points represent the 
path followed in the interaction zone during a semi-central 132 Sn+ 132 Sn 
collision at 1 AGeV (circles) and at 300 AMeV (crosses). 

In Fig. [3] we show the crossing density p cr separating nuclear matter from 
the quark-nucleon mixed phase, as a function of the proton fraction Z/A. We 
can see the effect of the (5-coupling towards an earlier crossing due to the larger 
symmetry repulsion at high baryon densities. In the same figure we report the 
paths in the (p, Z/A) plane followed in the cm. region during the collision of 
the n-rich 132 Sn+ 132 Sn system, at different energies. At 300 AMeV we are just 
reaching the border of the mixed phase, and we are well inside it at 1 AGeV . We 
expect a neutron trapping effect, supported by statistical fluctuations as well as 
by a symmetry energy difference in the two phases. In fact while in the hadron 
phase we have a large neutron potential repulsion (in particular in the NLpS 
case), in the quark phase we only have the much smaller kinetic contribution. 
Observables related to such neutron "trapping" could be an inversion in the 
trend of the formation of neutron rich fragments and/or of the ir~ /ir + , K°/K + 
yield ratios for reaction products coming from high density regions, i.e. with 
large transverse momenta. 

5 Perspectives 

We have shown that meson production in n-rich heavy ions collisions at in- 
termediate energies can bring new information on the isovector part of the 
in-medium interaction at high baryon densities. Important non-equilibrium ef- 
fects for particle production are stressed. Finally the possibility of observing 



precursor signals of the phase transition to a mixed hadron-quark matter at 
high baryon density is suggested. 
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